Harmonics generation and intense terahertz radiation from polar molecules at 
multiphoton resonant excitation in laser fields 
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The coherent radiation spectrum of two-level polar molecules with UV transition is studied at the 
multiphoton resonant excitation by a moderately strong laser field. The spectrum corresponding to 
harmonic generation and low-frequency radiation is investigated both analytically and numerically. 
Specifically, a mechanism for generation of intense smoothly tunable terahertz radiation by two- 
level molecular configurations (with dynamic Stark shifts) possessing permanent dipole moments, 
is considered. 

PACS numbers: 42.50.Hz, 42.65.Ky, 32.50.+d, 82.53.Kp 



I. INTRODUCTION 

The interaction of strong (moderately) laser radiation 
with the matter produces a number of significant pro- 
cesses that include the harmonic generation. The latter is 
of particular interest because of importance and potential 
applications of tunable, coherent, intense high-frequency 
radiation. Generation of high-order laser harmonics is 
a well-approbated process in the systems composed by 
atoms or small molecules and modeled most successfully 
by a " three-step recollision" physical picture (J, . The 
three-step theoretical models for high harmonic genera- 
tion does not include any participation of excited bound 
states in the emission of high harmonics and is essentially 
valid for harmonics with energy well above the ionization 
energy 

On the other hand, various studies have shown that 
for the certain systems such as quantum wells and large 
molecules [H, 0] , bound-bound transitions are more im- 
portant for harmonic generation than the coupling be- 
tween the ground and the continuum state [5|. This 
mechanism of harmonic generation without ionization 
can not provide high-frequency radiation, e.g. x-ray, 
but can be more efficient for generation of moderately 
high harmonics, in particular, for UV/VUV applications 
@. For description of strong light scattering process via 
bound-bound transitions, resonant interaction is of inter- 
est. Apart from its pure theoretical interest as a simple 
model, the resonant interaction regime may allow to in- 
crease considerably the efficiency of frequency conversion. 

Harmonic generation by a resonantly driven two-level 
atom has been studied in [7f . However, a two-level atomic 
system meets difficulties since the required for efficient 
multiphoton resonant excitation laser fields are so strong 
that many atomic levels and continuum will play impor- 
tant role [I]. So, two- level atomic system is not a good 
quantitative model for description of atomic response in 
case of multiphoton-resonant excitation. 

Nevertheless, in some atomic and especially molecular 
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systems one can avoid the difficulties with efficient multi- 
photon excitation based on the two-level model. As it has 
been shown the multiphoton resonant excitation 

of a quantum system subjected to a strong laser field is 
effective when the quantum system has permanent dipole 
moments (PDM) in the stationary states. Otherwise, the 
energies of the excited states of a three-level atom should 
be close enough to each other and the transition dipole 
moment between these states must be nonzero. Further- 
more, these systems have an advantage, which allow to 
generate radiation with frequency much lower than laser 
frequency [H|,[l6| ■ The effects of diagonal dipole elements 
on harmonic generation at the one-photon resonance has 
been considered in [l7|. However, for efficient genera- 
tion of moderately high harmonics by optical pulses one 
should consider systems with the bound-bound UV/VUV 
transitions and, therefore, the multiphoton resonant in- 
teraction regime is of special interest. The coherent scat- 
tering spectrum of a two-level system possessing PDM at 
multiphoton resonant excitation has been considered in 
Ref. [18|. However, the presented in this paper analytical 
results for radiation spectrum have only qualitative char- 
acter, in the meantime, the numerical calculations have 
been performed for such strong fields which break the 
condition of multiphoton resonant excitation (see below 
Eq. (H3])) and also could lead to detrimental ionization 
in real atomic/molecular systems. 

In the present paper, coherent light scattering by two- 
level polar molecules at multiphoton resonant excitation 
in the field of a moderately strong laser radiation is in- 
vestigated. It is shown that the study of this process 
is important for efficient population transfer, generation 
of moderately high harmonics, as well as for generation 
of low-frequency radiation, specifically, for realization of 
intense smoothly tunable terahertz radiation sources. A 
simple analytic expression for the time-dependent mean 
dipole moment, taking into account the dynamic Stark 
shifts, is obtained. In particular, based on this expres- 
sion results, concerning the main spectral characteristics 
of considering process, are in good agreement with results 
of performed numerical calculations. The effect of com- 
pensation of dynamic Stark shift on the radiation spec- 
trum is considered. We compare the spectra of two-level 
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systems with and without PDM and show the advantages 
of polar systems for efficient generation of desirable ra- 
diation. Special attention is paid to generation of coher- 
ent low-frequency radiation in such systems, particularly, 
smoothly tunable, intense terahertz radiation. 

It should be noted that the obtained results may be 
interesting for different type coherent light sources, such 
as polar gases, certain organic crystals or quantum dots 
and etc. [H, ll9l422T ]. Besides, all presented results can 
be scaled to other systems and diverse domains of the 
electromagnetic spectrum. 

The paper is organized as follows: in Sec. II we present 
the analytical model and derive the coherent contribution 
to the harmonic and low-frequency spectra. In Sec. Ill 
we present some results of numerical calculations of the 
considered issue without multiphoton resonant approxi- 
mation and compare the obtained spectra with analytical 
results. Finally, conclusions are given in Sec. IV. 



II. THEORY 

A. Basic Model 

We consider a two-level quantum system possessing 
PDM driven by a time-dependent, intense laser field. The 
Hamiltonian describing the interaction of the laser field 
with the considered quantum system is given within semi- 
classical dipole approximation by 

H = (e x + Vn) |1)(1| + (e 2 + V 22 ) |2)(2| 



+ (V 12 |l)<2|+h.c). 



(1) 



Here, e\ and £2 are the energies of the stationary states 
|1) and |2) of unperturbed molecule, £2 > E\- In Eq. (UJ 



Vqv = -d nu E COSLJ t, 



(2) 



is the interaction part of Hamiltonian with real ma- 
trix element of the electric dipole moment projection 
d nv = e • d,j„, and pump wave field is taken to be lin- 
early polarized, with unit polarization vector e, slowly 
varying amplitude Eq, and carrier frequency wo. The di- 
agonal terms in ([2|) describe the interaction due to the 
PDM and are crucial for effective multiphoton coupling. 

The technique which we employ in our calculations 
has been described earlier in Q. First we obtain dy- 
namic wave function of the system in the pump laser field, 
i.e., the solution of the time-dependent Schrodinger equa- 
tion with Hamiltonian (fT}. From these solution we de- 
duce formula for expectation value of the time-dependent 
dipole-moment operator. Note that dipole-moment form 
of radiation spectrum gives similar results with dipole- 
acceleration and dipole-velocity forms in case of long 
pulses [23j]. 

We consider Schrodinger equation: 



■ !*(*)) 
' at 



(3) 



with Hamiltonian (fT]) at multiphoton resonant excitation 
regime and assume that \S n \ <C luq, where 5 n is the n- 
photon resonance detuning, given by the relation: 



S n = 61 - e 2 + nuj . 



(4) 



Here and below, unless stated otherwise, we employ 
atomic units (ft = e = m e = 1). Our method of solv- 
ing © has been described in detail in [l(| and will not 
be repeated here. Under the generalized rotating wave 
approximation, the time-dependent wave function can be 
expanded as: 



|tf (t)) = e - l£lt {[3x(t) + ai(t)] |1> + [3a (t) + a 2 (t)] 



exp[-i(nujot + / (V22 - V n )dt)] 




|2) 



(5) 



where aj (i) are the time-averaged probability amplitudes 
and oti(t) are rapidly changing functions over pump wave 
period. The time-averaged amplitudes 3«(i) are: 



en = y^Cjj exp(iAji), 
3=1 



(6) 



where are the constants of integration determined by 
the initial conditions, and the factors Xj are the solutions 
of the second-order characteristic equation: 



A„ — A — F n 
—F n — A n — S n — A 



with the function: 



F n = -^-nu} J n (Z) 
dp 



and detuning: 



k 2 J 2 k (Z) 



(7) 



(8) 



(9) 



Here d p — d 2 2 ~ dn (let d p > 0) is the difference of 
PDM in two stationary states -in the excited and ground 
states. The argument of the ordinary Bessel function 
J n (Z) is the difference of dipole interaction energies in 
units of the pump wave photon energy: Z = d p Eo/u}Q. 
The terms F n and A„ describe the resonant coupling and 
dynamic Stark shift at n-photon resonance, respectively. 
In deriving these equations we have applied well-known 
expansion of exponent through Bessel functions with real 
arguments [Hj]: 



iZ sin a 
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J s (Z)e lsa . 



(10) 



Assuming smooth turn-on of the pump wave, the 
relation between the rapidly and slowly oscillating parts 
of the probability amplitudes can be written as: 



«i(t) = a 2 (t)— ^ -— , 



(11) 
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a 2 (i) = -ai(t)— > . (12) 



k — n 



It should be noted that the nonperturbative resonant ap- 
proach used here put the following restrictions: 



|F„|,|A„|,|£ n | <w 



(13) 



on the characteristic parameters of the considered prob- 
lem (for more details see [loll). As it was mentioned in 
Introduction, this condition is not fulfilled in Ref. (l8| . 



B. Harmonic generation 

We restrict our study to the coherent part of the spec- 
trum, i.e., the spectrum of the mean dipole moment 
which is predominant for the forward direction when the 
number of scatterers is relatively large [25|, [26j]. In the 
Schrodinger picture the coherent part of the spectrum is 
251: 



die 



W)) 



where 



{£>(*)> = <*(*)! §-a(o) m)) 



(14) 



(15) 



is the time-dependent mean dipole moment, which is the 
main observable quantity. With the help of the wave 
function ([5]) the expectation value of the dipole operator 
(TTSl) can be written as: 



<£>(*)> =dii|ai(t)| 2 +d 22 |a 2 (i)| 2 



^d 12 a\a 2 Jk (Z) e^-^ + c.c.j 



(16) 



where ai. 2 (i) = ai ^(t) + ai j2 (i). Combining Eqs. ([6]), 
(fTTjl , (fT2|) , and (|T6f one can calculate analytically the ex- 
pectation value of the dipole operator for an arbitrary ini- 
tial atomic state. The Fourier transform of (D(t)) gives 
the coherent part of the dipole spectrum. 

The solution © for the system initially situated in the 
ground state, is: 



(71 
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which expresses the Rabi oscillations with frequency 



(17) 
(18) 



n n = + (2A 



The generalized Rabi frequency at n-photon resonance 
has nonlinear dependence on the amplitude of a pump 
wave field. As we can see from (fT7|l . ([T8]) the amplitudes 
of Rabi oscillations are diminished due to dynamic Stark 
effect. The destructive effect of the latter can be com- 
pensated by an appropriate detuning: 



-2A, 



(20) 



Replacing the probability amplitudes in (|16[) by the cor- 
responding expressions (flT)) and (jl8]l . one can derive 
the final analytical expression for (D(t)). As is seen 
from Eqs. flSJ, ©, the multiphoton coupling is pro- 
portional to the ratio dyz/d p , while the dynamic Stark 
shift is proportional to d^/dp. Since large dynamic 
Stark shifts make difficult the maintenance of consider- 
able population transfer, here we consider systems with 
A* = I di2 1 /dp « 1. Taking into account the smallness of 
the parameter /i, from (fTBf in the first order of approxi- 
mation by the small parameter /i we obtain the following 
compact analytic formula: 

(D(t)) = Co + J2[S k sin(fcwoi) + C k cos(ku t)}, (21) 



where 

C = d n + 



, 2^ 



2F n 2A r 



-Jn{Z) 



X (1 - cos(Q„*)), (22) 
describes the low-frequency part of the spectrum, and 
2F n nJ n+k (Z) 



S k = d 



12 



■ sin(fi„t), 



(23) 



2F n 2A» + S n nJ n+k (Z) 



(1 - cos(fi»t)), (24) 



describe harmonic radiation. As we can see from (|23p~ 
(|24p . the harmonic spectrum consists of triplets, with 
harmonic and two hyper-Raman lines with frequencies 
displaced by ft n . From (fl9|) follows that frequencies of 
hyper-Raman and low-frequency lines can be smoothly 
tuned via variation of the laser field strength and detun- 
ing. 

The expression for (D(t)) shows that intensities of the 
harmonics mainly determined by the behaviour of Besscl 
function. Since the Bessel function J m (Z) steeply de- 
creases with increase of index m > Z , then the cut-off 
harmonic s c is determined from the condition s c — n ~ Z. 
Then, from the estimation for cut-off harmonic follows 
that upper limit of the frequency which can be effec- 
tively generated by direct n-photon excitation, is higher 
for the systems with larger difference of energy and dipole 
moments in stationary states (uj c ~ e 2 — E\ + dpEg) . Al- 
(19) though cut-off frequency is not directly depend on the 
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FIG. 1: (Color online) The dependence of the amplitude of 
state-populations oscillation A — 4F^/f2^ on parameter Z 
for zero detuning (S n = 0) and d\2jd v = 0.1. The solid (red) 
line corresponds to the three-photon (n = 3) resonance; the 
dashed (green) line corresponds to the five-photon (n — 5) 
resonance; the dotted (blue) line corresponds to the eight- 
photon (n = 8) resonance. 



hyper-Raman lines becomes weaker than the other two 
lines of the triplet. Distribution of intensities is essen- 
tially changed for almost complete population transfer, 
i.e., when f2„ m 2\F n \ ^> |2A„ + <5 n |. In this case the 
intensities of the hyper-Raman components become on 
the same order and harmonic component practically van- 
ishes. The latter is proportional to /i 2 and appears only 
in the next order of approximation by the small param- 
eter jJL. 



C. LOW-FREQUENCY RADIATION 

As it has been mentioned above, the expectation value 
of the dipole operator (|2Tj) also has low-frequency term 
{D\ w{t)) which is given by the formula: 



(Aow(*)) = - d P 



IFl d 12 2F n 
fi 2 HQ 



2 At, 



a, 



-Jn (Z) 



COS (fi„t) , 



(26) 



photon order n of the resonance, but taking into account 
harmful ionization, the large n is preferred for efficient 
harmonic generation at the given strength of the laser 
field. Note that due to the oscillating character of Bessel 
function, harmonics can be strengthened or weakened via 
variation of the laser field strength within the resonance 
width. 

The intensity of the harmonic components is propor- 
tional to F^/fl n , which is the amplitude of population 
oscillations of the states. Thus, for the effective radia- 
tion generation one should provide a considerable popula- 
tion transfer between the ground and excited states, i.e., 
|F„| ~ Q n . As it follows from ©, (|9|), and ([19]), besides 
the ratio d\ijd v ^ parameter Z may also have substantial 
effect on population transfer. In Fig. 1 we plot depen- 
dence of the amplitude of population oscillations of the 
states on parameter Z , for zero detuning and different 
values of the photon order n. The complete population 
transfer is achieved when A„ = 0. In fact, for consid- 
erable population transfer it is important to compensate 
dynamic Stark shift by an appropriate detuning, specifi- 
cally for high-order multiphoton excitation in the weaker 
laser fields. 

The pattern of dipole intensities in harmonic triplet is 
determined from the expressions (|2U1) - (|2"4"|) : 



Jn-k (Z) + J n +k (Z) Q n 
Jn-k (Z) — J n +k (Z) 2A„ + 5„ 



-i 2 



(25) 

In general, due to oscillating character of the Bessel func- 
tions, the triplet pattern can be significantly changed de- 
pending on the pump wave field strength. In moderately 
strong laser fields Z < 1 and in case of considerable pop- 
ulation transfer, i.e., when Sl„ ~ |2A„ + <5„|, one of the 



according to Eq. (|22j) . Formula (|26|) describes the ra- 
diation at the smoothly tunable frequency Q n , i.e., at 
the frequency of oscillations of the population inversion. 
This low-frequency radiation also depends on the ampli- 
tude of population oscillations of the states. According 
to Eq. (|2"6"|) . with increase of the amplitude of Rabi oscil- 
lations the intensity of low-frequency radiation increases. 
When population transfer is almost complete, e.g., when 
dynamic Stark shift is compensated by appropriate de- 
tuning, the low-frequency term can be written as: 



(AowW) - — ^COS(fiflt), 



(27) 



where the Rabi frequency Qr is determined by the ex- 
pression: 



0/ 



2-^nujQ J n (Z) 
d p 



(28) 



In this case the emitted radiation intensity does not de- 
pend on the pump wave field strength and is determined 
only by the difference of PDM in the excited and ground 
states. 

For the considered problem, the Rabi frequency may be 
varied within the following limits: 27r /T r << CIr << luq, 
where T r describes the relaxation time in the molecular 
system. Hence, for electronic transitions FIr can lie from 
microwave to infrared frequencies. Here, we concentrate 
on terahertz region of electromagnetic spectrum (0.1 — 10 
THz) due to the important properties of such radiation 
and its numerous applications for sjjectroscopy, sensing, 
imaging, and etc (for a review see [271] ). 

For comparison of considered scheme with the other 
schemes of THz radiation, we make some estimations for 
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total radiation power of the ensemble of N molecular 
emitters. It is assumed that the molecules are oriented 
in the same direction. Though, in case of freely rotating 
molecules temperature-dependent distribution of molecu- 
lar orientations reduces transition and permanent dipole 
moments, the intensity of emitted radiation also will be 
reduced to some extent. Due to phase-matching factor 
the coherent radiation entirely occurs almost in forward 
direction. The power of low- frequency radiation was eval- 
uated for resonantly driven three-level atomic-molecular 
ensemble in [l6[ and will not be repeated here in detail. 
We will adopt the estimations made for circular cylinder 
model of superradiant emitter (in a limit of long needle) 
for ensemble of 2-level polar molecules. The radiation 
power for ensemble of polar molecules coherently excited 
by Gaussian laser beam with a waist wq can be approxi- 
mated by formula: 



16c 2 



(29) 



where is a degree of monochromaticity specified by 
driven field variation in the interaction region, L c is the 
coherent length over which the low-frequency can be built 
up coherently and Nq is the density of polar emitters. 
The corresponding incident pulse duration should be: 
t ^3 27rn/(f2fl Aji). As we can see from (|29[) . it is more de- 
sirable to generate THz radiation via low-order multipho- 
ton excitation. Here we make estimations for two-photon 
resonance (n — 2); the incident laser beam waist is taken 
to be: wq 1 mm, the coherence length: L c ^3 10 cm, 
the difference of PDM in the excited and ground states: 
d p ~ 5 a.u. and emitters density: No ~ 10 15 cm -3 . For 
the moderate monochromaticity An = Silu/ftu ~ 0.1 at 
SIr/(2tv) ~ 3 THz, the incident pulse duration should be: 
T 6 ps. Estimated total radiation power P ~ 100 W 
is comparable with parameters for most powerful THz 
sources. Note that submillimeter-sized array of reso- 
nantly driven 2-level quantum dots with induced dipole 
moment ~ 10 D gives the power of microwatt level |15| . 

In case of a solid medium, the radiation power could 
be larger due to higher densities of matter. The other ad- 
vantage of the solid medium is the degree of anisotropy, 
since molecules can practically be oriented in the same 
direction, e.g., in some organic crystals or when the 
molecules are inserted in a solid matrix that is transpar- 
ent to the radiation. The volume of coherent radiation 
for a solid medium is ~ A 3 , and the number of polar emit- 
ters can be estimated by N ~ A 3 No, which at densities 
iVo ~ 10 21 cm -3 gives N ~ 10 15 . The rough estimations 
show that using submillimeter-sized crystals with ~1% 
ordered dipole moments of emitters one can achieve gi- 
ant total radiation power for coherent THz radiation up 
to megawatt level. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section we present numerical solutions of time- 
dependent Schrodinger equation for a two-level model 
with Hamiltonian (JTJ. The set of equations for the 
probability amplitudes has been solved using a standard 
fourth-order Runge-Kutta algorithm [28j . The Fourier 
transformations for estimation of power spectra are per- 
formed using the fast Fourier transform technique. For 
smooth turn-on of pump wave field, the latter is de- 
scribed by the envelope in hyperbolic tangent tanh(t/r) 
form, where r characterizes the turn-on time and cho- 
sen to be 207r/a;o. The transition frequency is set to be 
£2 — £1 = 0.2 a.u. and lies in UV region, but the results 
can be scaled to different energy regimes. The transition 
dipole moment is chosen to be d\i = 0.5 a.u., and differ- 
ence of PDM in the excited and ground states: d p = 5 
(except where it is declared otherwise). It is as- 



a.u. 

sumed that initially all the population is in the ground 
state |1) . 

Figure 2 shows dipole spectrum (coherent part) as a 
function of harmonic order for four-photon (n — 4) reso- 
nant excitation. The pump held strength is E = 0.01 a.u. 
and multiphoton detuning 5 n is set to zero. The solid 
(red) line corresponds to numerical calculations, while 
the dashed (green) line corresponds to the approximate 
expression (1211) . For better visibility, the spectrum cor- 
responding to analytical calculations has been slightly 
shifted to the right. As we can see from the Fig. 2, the 
analytical formula (|21[) is in good agreement with the 
exact results. For employed field strengths the triplet 
structure of harmonics and low-frequency line are not 
clearly seen, hence the latter is illustrated in the inset. 

Dipole spectrum as a function of harmonic order for 
five-photon (n — 5) resonant excitation with almost com- 
plete population transfer is plotted in Fig. 3. The pump 
field strength is E = 0.01 a.u. and detuning for com- 
pensation of dynamic Stark shift is 5 n = 0.00014 a.u.. 
As in Fig 2, the solid (red) line corresponds to numerical 
calculations and the dashed (green) line corresponds to 
the approximate expression (1211) . The numerical results 
confirm the analytical expression for expectation value 
of the mean dipole moment. These figures confirm the 
estimation s c — n ~ Z for cut-off position. 

In Fig. 4 we plot characteristic harmonic triplet (here, 
for 6th harmonic) under the same conditions of excita- 
tion, as in Fig. 3. The solid (red) line corresponds to de- 
tuning, compensating dynamic Stark shift, and dashed 
(green) line corresponds to zero detuning. For better 
visibility, the spectrum corresponding to compensating 
detuning has been slightly shifted to the right. The nu- 
merical results confirm estimations for lines position in 
harmonic triplet and the fact that via compensation of 
dynamic Stark shifts one can attain the higher intensities. 
Distribution of intensities, particularly, damping of har- 
monic lines in case of compensating detuning and lower- 
frequency hyper-Raman lines at zero detuning, have been 
predicted in section II. 
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FIG. 2: (Color online) The logarithm of the coherent part 
of the spectrum Sc(w) at four-photon (n = 4) resonance 
with zero detuning (8 n = 0). The energy level difference 
is £2 — £i = 0.2 a.u., laser field strength E — 0.01 a.u., and 
dipole moments are di2 = 0.5 a.u. and d p = 5 a.u.. The solid 
(red) line corresponds to numerical calculations; the dashed 
(green) line corresponds to the approximate solution (for bet- 
ter visibility the latter has been slightly shifted to the right). 
The inset shows the low-frequency peak at Rabi frequency. 
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FIG. 3: (Color online) The logarithm of the coherent part 
of the spectrum Sc(^) at five-photon (n = 5) resonance. 
The energy level difference is £2 — £i = 0.2 a.u., detuning 
5 n — 0.00014 a.u., laser field strength E = 0.01 a.u., and 
dipole moments are d\2 = 0.5 a.u. and d p = 5 a.u.. The solid 
(red) line corresponds to numerical calculations; the dashed 
(green) line corresponds to the approximate solution (for bet- 
ter visibility the latter has been slightly shifted to the right). 



To demonstrate the effects of PDM we plot the spec- 
trum of the two-level system without PDM under the 
same conditions of excitation, as in Fig. 3. As is shown 
in Fig. 5, there are fewer harmonics and with negligibly 
small amplitudes. Practically the full radiation is con- 
centrated on the incident radiation frequency. Even har- 
monics, as well as the low-frequency radiation are absent 
because of inversion symmetry. Note that these essen- 



FIG. 4: (Color online) The structure of satellite peaks around 
the 6th harmonics under the same conditions of excitation as 
in Fig. 2. The solid (red) line corresponds to zero detuning 
(8 n = 0); the dashed (green) line corresponds to detuning 
S n = 0.00014 a.u. (for better visibility the latter has been 
slightly shifted to the right). 




4 6 
Harmonic order 

FIG. 5: Same as in figure 2, but for the case of zero dipole 
moments in stationary states. 



tial distinctions in intensities of the harmonics are for 
the pump field strengths with negligible ionization and 
dynamic Stark shifts. To achieve the efficient harmonic 
generation when d p = 0, required laser fields should be 
comparable to characteristic fields of considered system: 
Eq > (£2 ~ £1) /d\2, which will cause complete ioniza- 
tion. 

We also made numerical calculations for low-frequency 
part of the radiation spectrum. Figure 6 displays the 
dependence of the emission rate on emission frequency 
at five-photon resonance. The lower line corresponds to 
zero detuning and upper line to detuning compensating 
dynamic Stark shift. The calculations have been made 
for the same field strengths. As we can see from the Fig. 
6, the frequency and intensity of the emission lines ap- 
preciably depend on the magnitude of the detuning and 
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FIG. 6: The low- frequency emission rate as a function of the 
emission frequency under the same conditions of excitation as 
in Fig. 2. (a) upper line corresponds to compensated dynamic 
Stark shift, (b) lower line corresponds to zero detuning [5 n = 
0). 



compensation of the Stark shift is of considerable impor- 
tance. The numerical results for spectral characteristics 
are in good agreement with estimations based on (|26[) 
and ((27) . 

It should be noted that evenly charged molecular ions 
at large internuclear distances lilt , driven flux qubits 
(29| , and hydrogenlike atoms |l0l.ll2l [l3| are good exam- 
ples of efficient direct multiphoton excitation. The men- 
tioned three-level configurations are also good candidates 
for moderately high harmonics generation [6(. However, 
in contrast to 2-level model with PDM, due to inver- 
sion symmetry the spectrum of three-level system does 
not contain even harmonics (with corresponding hyper- 
Raman lines) and low-frequency components under the 
conditions of excitation considered here. 



IV. CONCLUSION 



direct multiphoton resonant excitation. With the help 
of nonperturbative resonant approach, we have studied 
approximate stationary solution of the Schrodinger equa- 
tion and obtain a simple analytical expression for scat- 
tered coherent radiation. We have also numerically inves- 
tigated the harmonic spectra and effects of compensation 
of the dynamic Stark shift on emission spectrum. As has 
been shown, the compensation of dynamic Stark shift 
by appropriate detuning could substantially improve ef- 
ficiency of radiation generation. The numerical results 
are in good agreement with obtained analytical results 
for expectation value of dipole operator. We have also 
shown that the existence of PDM enables the consider- 
able population transfer in moderately strong laser fields, 
and therefore could contribute to generation of higher 
harmonics with sufficiently large amplitudes. In addi- 
tion, due to broken inversion symmetry of systems with 
PDM there is substantial emission at frequency much 
smaller than the laser frequency. The considered scheme 
may serve as a promising method for efficient generation 
of moderately high harmonics and intense widely tunable 
low-frequency radiation. Note that for coherent radiation 
of an molecular ensemble one needs more rigorous treat- 
ment, which should account for degree of anisotropy in 
the orientational distribution of polar molecules, as well 
as collective effects in the medium. Work in this direc- 
tion is in progress and will be presented in forthcoming 
paper. 
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